Baryon phase-space density in heavy-ion collisions 
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The baryon phase-space density at mid-rapidity from central heavy-ion collisions is estimated from 
proton spectra with interferometry and deuteron coalescence measurements. It is found that the 
mid-rapidity phase-space density of baryons is significantly lower at the SPS than the AGS, while 
those of total particles (pion -I- baryon) are comparable. Thermal and chemical equilibrium model 
calculations tend to over-estimate the phase-space densities at both energies. 
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Baryon phase-space density affects the dynamical evo- 
lution of heavy-ion collisions from the initial stage to final 
freeze-out. It is one of the essential ingredients of ther- 
mal and chemical equilibrium models often used in study- 
ing heavy-ion collisions |jl|,|| . The studies of baryon and 
meson phase-space densities may provide information on 
particle freeze-out conditions and entropy production in 
heavy- ion collisions Experimentally, phase-space 

densities can be estimated from particle yields and source 
sizes inferred from two-particle interferometry (^-Ml and, 
in the case of nuclcons, from deuteron coalescence 
The values of the phase-space densities may be then com- 
pared to Bose- Einstein or Fermi-Dirac statistics to assess 
the degree to which thermal and chemical equilibria are 
established [9[-[Tl|. One unique advantage of such com- 
parisons is that they do not require the extrapolation of 
data to all of phase-space . 

The freeze-out particle spatial density is also of inter- 
est. However, it is known that, in high energy heavy- 
ion collisions, the freeze-out space point of a particle is 
correlated with its momentum. Because of this space- 
momentum correlation, the size parameters extracted 
from two-particle interferometry at low relative momen- 
tum do not reflect the full size of the particle-emitting 
source jlj]. Hence, one cannot readily deduce the parti- 
cle spatial density from the size parameters and particle 
yields which is from the entire source. 

There have been several studies of the freeze-out phase- 
space density of pions flpl-pdl. In this paper, we study 
the freeze-out phase-space density of baryons. We first 
estimate the baryon phase-space density at mid-rapidity 
in heavy ion collisions at the AGS and SPS using two- 
proton interferometry data. We then investigate the ratio 
of deuteron to proton invariant cross-sections as a means 
to extract the baryon phase-space density Finally, 
we compare our estimates to thermal and chemical equi- 
librium model calculations. 

We start from the definition of particle phase-space 
density, 

/P,x) = -— -, (1 
dpdx 

which is Lorentz invariant. Consider identical particles 
in a small momentum cell, (p,p -I- dp), which occupy 



a spatial volume V^(p) in the particle rest frame. (We 
specify V{p) in the particle rest frame because available 
two-proton correlation measurements p^ , p^ are done in 
such frame.) If these particles are uniformly distributed, 
then the spatial averaged phase-space density is 



(/(P)> = 



E dN/dp 
m V{p) 



(2) 



where m is the particle rest mass, and E/m is due to the 
Lorentz boost. 

Instead of a uniform distribution in space, if we take 
a Gaussian form for the particle density profile (as often 
used to extract source size parameters from two-particle 
interferometry 1 17 jl^ ) , 



/(P,x) 



E dN/dp 
^ [V2^RGip)y 



[x-xo(p)]V2fl.?;(p) 



then Eq. (0) becomes 



(/(P)> = 



/ /2(p,x)dx _ E dN/dp 



//(p,x)dx m[20FiJc(p) 



(3) 
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Rg is a measure of the relative separation of two particles 
close in momentum in the pair rest frame when the later 
particle freezes out |l9| . In the above equations, however, 
it has been assumed that the particles freeze-out at an 
instant time, and the phase-space density is calculated at 
that time. When the freeze-out process has a finite time 
duration, the phase-space density should be considered 
as a "time-averaged" quantity. 

Without knowing the exact form of the space- 
momentum correlation, xo(p), one cannot perform the 
integration over p. Hence, one cannot obtain the spatial 
density. However, the operation can be carried out over 
the spatial coordinate x leading to Eq. (^) . In the follow- 
ing discussions, we will use the Gaussian density profile 
and denote V{p) = [2^RGip)]^ ■ 

The averaged phase-space density weighted by particle 
number density is given as 



r 1 f T^dN\'^ dp 

j{f{p))dN _J T^n^[^ii^) -E 



JdN 



(5) 
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We are mainly interested in the baryon phase-space den- 
sity. Since the transverse momentum (px) dependence 
of Ro for mid-rapidity proton source has not been mea- 
sured, we assume an pT independent volume V{y,pT) = 
V{y) in the following discussions. Taking an exponen- 
tial function in iriT = y'm^ + p^ for the invariant cross- 



section 



ifiy)) 



dN/dy 
2TTT(T+m) 

dN/dy 



^-(mT-m)/T obtain 



oo ~2{mT-rn)/T 



27rmT2(T + m)2 7„ V{y,pT) 

T + 2m dN/dy 
87rmT{T + my V{y) 



-mTdrriT 
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Bertsch obtained similar results for mesons from two- 
particle correlation formula , assuming an exponential 
transverse distribution in pT. Our results above are valid 
for both mesons and baryons, and are consistent with 
Bertsch's if same assumptions are made. 

Using Eq. (||) and available data |l^ , |20| -|26|| , the av- 
erage phase-space densities are estimated for freeze-out 
pions, protons and baryons at mid-rapidity in central 
Au+Au or Pb+Pb coUisions at the AGS and SPS. The re- 
sults are listed in Table |. The statistical and systematic 
errors on the phase-space densities are estimated from 
the corresponding errors on dN/dy, T, and Rq. In es- 
timating the baryon phase-space density, it is assumed 
that the source sizes of other baryons are the same as 
proton's. 

The source size parameters Rg are estimated from 
pion and proton interferometry measurements in invari- 
ant relative momentum p^pl| , p5|,p6|| . For pions, Rg{p) 
is observed to decrease with pT [p5|-p8[ due to the space- 
momentum correlation. However, we have only used the 
low pt pion source size in our estimation. This intro- 
duces an underestimate of the pion phase-space density, 
for instance, by about 25% for the SPS data. For the 
mid-rapidity proton source size at the AGS, we have used 
the two-proton correlation measurement at forward ra- 
pidity jlj] and estimated as Rq — 3.5 ± 0.5 (syst.) fm 
such that the low end of the systematic error equals to 
the initial nucleus size Rq — 3.0 fm. This estimate is 
consistent with a hypothetical l/y/mr scaling from pion 
interferometry measurements. The errors introduced by 
assuming a pt independent volume are not included in 
the systematic errors quoted in Table |. 

It is known that the pion and proton transverse dis- 
tributions cannot be satisfactorily described by a single 
rriT exponential p0|p2| -|2^. In order to check the errors, 
we used the available data to calculate the phase-space 
density by Eq. (||) and the results are consistent with 
the above estimates within 15%. This is also reflected 
in the systematic errors on the inverse slope parame- 
ter (T) quoted in Table |. In fact, the values of the 
phase-space density extracted from this work are consis- 
tent with those obtained in [p|-pT[|. 



As seen in Table |, the freeze-out pion phase-space 
densities are similar (~0.1) for both AGS (-y/s « 5 
AGeV) and SPS (^/i w 20 AGeV) energies. But the 
proton(baryon) phase-space density is significantly lower 
than the pion's at both beam energies. On top of that, 
the mid-rapidity baryon phase-space density decreases as 
the beam energy increases. The similarity of the total 
phase-space density (pion -|- baryon) at freeze-out per- 
haps indicates that, in high energy collisions, the freeze- 
out process occurs at a fixed condition (local tempera- 
ture and density) which is determined by the free cross 
sections as argued by Pomeranchuk in 1951 |^9|. Teva- 
tron data seem to indicate a constant particle density 
in pp annihilations jsoj. Note that at the AGS energy, 
pion to proton ratio is about one so the baryon density 
is the relevant control parameter for freeze-out At 
the SPS energy, the pion to proton ratio is about six 
making mesons the dominant player at freeze-out. This 
will be certainly true in the upcoming RHIC collisions at 
Vs = 200 AGcV. 

We have assumed constant ^(p) in estimating (/) 
above. This is because, as mentioned before, a px depen- 
dent two-proton correlation measurement is not presently 
available in high energy nucleus-nucleus collisions due to 
the lack of statistics. Such measurement, however, may 
be obtained from deuteron coalescence [^2|-[3^]. 

In the coalescence model with a Gaussian density 

profile, RUp,) = I-'" /{^M 
where the factor 3/4 comes from spin consideration, and 
subscripts 'p' and 'd' denote proton and deuteron, re- 
spectively. From Eq. (0), we obtain 



(/(Pp)> 



4 [Ed dNi 



3(27r)3 V"id dpd/p^^app/ \"^p ^Pp 



^p dN^ 



(7) 



where the neutron and proton differential cross sections 
are assumed to be identical []. The average proton phase- 
space density [Eq. (||)] can be approximated as (/) = 
g^5^ / dN^ j J dNp . Therefore, the pT averaged proton 
phase-space density is 



ifiy)) = 



1 



{dN/dy)d 



6(2^)3 {dN/dy)p' 
and that averaged over the whole phase-space is 

1 iVd 



6(27r)3 Np 



(8) 



(9) 



^ Various models | p9|j40| show that the neutron to proton 
dN/dy ratio at mid-rapidity in Au-|-Au (AGS) and Pb-|-Pb 
(SPS) central collisions is about 1.05. 
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as pointed out long ago by several authors 

Using Eq. (||) and available deuteron data | |21| , [41| -[43[ 
we estimate the proton phase-space density at the AGS 
and SPS to be (/)(27r)3 = 0.015 ± 0.001 (stat.) and 
0.0029 ± 0.0004 (stat.), respectively. These results are in 
a good agreement with those obtained using two-proton 
interferometry data (Table |). 

Now we compare our results to calculations by local 
thermal and chemical equilibrium models. In these mod- 
els, temperature (Tq) and baryon chemical potential (/is) 
at chemical freeze-out are extracted from particle ratios 
[p],pp^-p^ ; temperature at kinetic freeze-out is extracted 
from particle spectra with certain assumptions about col- 
lective flow 1 25 3^ . These results as obtained from AGS 
and SPS data are reproduced in Table ||. 

In local thermal and chemical equilibrium models, the 
phase-space density is {f{E)) = (27r)"3 / [e^/^" - l] for 
pion and {f{E)) = 2(27r)-3 / [e(£;-A'B)/To + i] for pho- 
ton, respectively. Here E is the particle energy, Tq the 
temperature, and fj,B the baryon chemical potential. The 
average particle phase-space densities at mid-rapidity, 
(/) = ^{f{E))'^dpl / ^{f{E))dpl, are calculated sepa- 
rately for chemical and kinetic freeze-out. The results are 
listed in Table 0. The phase -space densities at chemical 
freeze-out are higher by a factor of 2.7 ± 1.2 (AGS) and 
4.0 ±1.3 (SPS) than estimated from experimental data 
(Table H). Likewise, the phase-space densities calculated 
for kinetic freeze-out are higher by a factor of 1.5 ± 0.7 
(AGS) and 2.2 ± 1.2 (SPS). 

One important assumption in local thermal and chem- 
ical equilibrium models is that all particles freeze-out at 
the same time, and they occupy the same volume so that 
the volume is divided out in particle ratios. From this 
assumption and the fact that particle numbers stay con- 
stant after chemical freeze-out, an extrapolation has been 
made for the chemical potential from chemical to kinetic 
freeze-out. This is equivalent to the assumption that 
the system undergoes isentropic expansion after chem- 
ical freeze-out [^,0 [[ After chemical freeze-out, the 
system expands and the temperature decreases. This 
leads to a factor of 3 increase in the volume at both the 
AGS and SPS energies. Taking the average radial flow 
velocity as Vi « 0.5c one may estimate the corre- 
sponding expansion time: (v^— l)RA/vi ~ 6 fm/c with 
Ra = 1.2^1/3. 

In summary, we have studied the mid-rapidity baryon 
phase-space density at freeze-out through (1) proton 
yield and two-proton interferometry measurements, and 
(2) the ratio of deuteron to proton differential cross- 



* In this procedure, it is explicitly assumed that local ther- 
mal and chemical equilibrium is always maintained for pions 
and protons during expansion of the collision system. These 
pions include resonance decay products. 



sections. For heavy-ion collisions at both the AGS and 
SPS, the estimated phase-space density of baryons is sig- 
nificantly lower than of pions. The baryon phase-space 
density at the AGS is about three times that at the SPS, 
while the values for the pion phase-space density are sim- 
ilar at both energies. The thermal and chemical equi- 
librium model calculations are higher than experimental 
data at both beam energies. 

Note that the present analyses suffer from large sys- 
tematic uncertainties. In order to obtain more precise 
results for baryons, the proton interferometry or deuteron 
coalescence data as a function of the transverse momen- 
tum are needed. In addition, we only studied the particle 
phase-space density at mid-rapidity. It will be interest- 
ing to repeat the study at projectile/target rapidities and 
compare with thermal model predictions. 
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TABLE I. The average phase-space density {/) estimated for freeze-out pions, protons and baryons at mid-rapidity in 
central Au+Au or Pb+Pb collisions at the AGS (upper section) and SPS (lower section). Input parameters are mid-rapidity 
dN/dy, inverse slope T, and Gaussian source size Rg in the pair rest frame. Unless specified, for each quantity the first error 
is statistical and the second error is systematic. 



particle dN/dy 


r (MeV) 


Rg (fm) 


(/)(2^)' 


7r+ f» TT" 62 ± 2 ± 8 [E 
proton 61 ± 1 ± 8 |o 
baryon 140 ± 2 ± 20 | 


170 ± 2 ± 10 1 
270 ± 3 ± 10 i 


20[ 6.2 ±0.4 ±0.4 (21 
20| 3.5 ± 0.5 (syst.) 
// 


0.087 ±0.017 ±0.022 
0.014 ± 0.006 (syst.) 
0.032 ± 0.015 (syst.) 


n+^TT- 180 ±10 ±20 
proton 27±4±5 
baryon 68 ± 8 ± 10 


[22 
-~ 

^ E 


|23| 206 ± 2 ± 20 1 

1 285 ± 5 ± 30 i 
4) // 


22 
22 


M 7.4 ± 0.4 (stat.±syst.) 1^,^ 0.106 ± 0.028 (stat.±syst.) 

^ 3.85±0.15i;5:25 111 0.0043 ± o.ooo8i;5 [];523 

// ' 0.0109 ± 0.0018+o'oo5g 



TABLE IL The average mid-rapidity phase-space densities for pions and protons at chemical ((/)ch) and kinetic ({/)kin) 
freeze-out from the thermal and chemical equilibrium model. Input parameters are temperatures (Tq in MeV) at chemical and 
kinetic freeze-out, and baryon chemical potential {fiB in MeV) at chemical freeze-out extracted from central Au+Au or Pb±Pb 
experimental data at the AGS (upper section) and SPS (lower section). All errors are statistical only. See text for estimation 
of at kinetic freeze-out. 



Chemical freeze-out 



Kinetic freeze-out 
92 ±2 



particle 



(/>ch(2^) 



(/>ki„(2^)^^ 



To = 125 
MS = 540 ± 7 




To 



590 ± 5 



IT ~ TV 

proton 



n 1 79+U.()U« 



0.107 ±0.004 
0.021 ±0.002 



To = 168 ±2.4 Ml 
/is = 266 ± 5 |45[| 



To = 120 ± 12 |25| 

fiB = 389 ±4 



TT « TT 

proton 



0.249 ±0.004 
0.017 ±0.001 



0.162 ±0.023 
0.0096 ± 0.0043 
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